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ABSTRACT 

Context. Hydrogen-rich, DA-type white dwarfs are particularly suited as primary standard stars for flux calibration. State-of-the-art 
NLTE models consider opacities of species up to trans-iron elements and provide reliable synthetic stellar-atmosphere spectra to 
compare with observations. 

Aims. We will establish a database of theoretical spectra of stellar flux standards that are easily accessible via a web interface. 
Methods. In the framework of the Virtual Observatory, the German Astrophysical Virtual Observatory developed the registered ser- 
vice TheoSSA. It provides easy access to stellar spectral energy distributions (SEDs) and is intended to ingest SEDs calculated by any 
model-atmosphere code. In case of the DA white dwarf G191-B2B, we demonstrate that the model reproduces not only its overall 
continuum shape but also the numerous metal lines exhibited in its ultraviolet spectrum. 

Results. TheoSSA is in operation and contains presently a variety of SEDs for DA-type white dwarfs. It will be extended in the 
near future and can host SEDs of all primary and secondary flux standards. The spectral analysis of G191-B2B has shown that our 
hydrostatic models reproduce the observations best at T s g = 60 000 ± 2000 K and log g = 7.60 ± 0.05. We newly identified Fe VI, Ni VI, 
and ZnlV lines. For the first time, we determined the photospheric zinc abundance with a logarithmic mass fraction of -4.89 (7.5 x 
solar). The abundances of He (upper limit), C, N, O, Al, Si, O, P, S, Fe, Ni, Ge, and Sn were precisely determined. Upper abundance 
limits of about 10% solar were derived for Ti, Cr, Mn, and Co. 

Conclusions. The TheoSSA database of theoretical SEDs of stellar flux standards guarantees that the flux calibration of all astronom- 
ical data and cross-calibration between different instruments can be based on the same models and SEDs calculated with different 
model-atmosphere codes and are easy to compare. 

Key words, standards - stars: abundances - stars: atmospheres - stars: individual: G191-B2B - virtual observatory tools - 
white dwarfs 


1. Introduction 

In the framework of the Virtual Observatory (VO), the German 
Astrophysical Virtual Observatory (GAVO) project provides 
synthetic stellar spectra on demand via the registered Theoretical 
Stellar Spectra Access (TheoSSA) VO service (Rauch 2008a; 
Rauch & Nickelt 2009; Rauch et al. 2009). These SEDs can be 
used for spectral analyses (Rauch et al. 2010; Ringat & Rauch 
2010; Rauch & Ringat 2011; Ringat et al. 2012) or serve as 


* Based on observations with the NASA/ESA Hubble Space 
Telescope, obtained at the Space Telescope Science Institute, which is 
operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26666. 

** Based on observations made with the NASA-CNES-CSA Far 
Ultraviolet Spectroscopic Explorer. 

*** Figures 1, 6, 10-12, 23, A. 1, A.2 and Tables 2^1 are available in 
electronic form at http : / /www . aanda . org 

**** Table 5 and Figs. A. 1 and A.2 (FITS files) are only available at the 
CDS via anonymous ftp to 
cdsarc .u-strasbg . fr (13® . 79 . 128 . 5) or via 
http ://cdsarc .u-strasbg . fr/viz-bin/qcat?J/A+A/56S/Al®6 


ionizing spectra for e.g. photoionization models of ionized neb- 
ulae. The registered TMAW VO tool 1 , that allows to calculate 
individual NLTE model atmospheres considering opacities of 
H. He, C, N, O, Ne, Na, and Mg, provides additional SEDs 
which are automatically ingested by TheoSSA. Figure 1 shows 
the complete action scheme for a VO user to retrieve an SED. 

With the increasing usage of TheoSSA over the last years, 
it became necessary to demonstrate the reliability of the SEDs. 
We established simple benchmark tests (Ringat et al. 2012) to 
show the achievable analysis precision, e.g. in the determina- 
tion of effective temperatures (T e a) and surface gravities (loggr), 
in cases that TMAW SEDs are used which are calculated with 
standard model atoms that are limited in the number of atomic 
levels treated in NLTE. This guarantees model calculations in a 
reasonable time for a VO user. 

Since 2012 TheoSSA also includes synthetic spectra of spec- 
trophotometric standard stars. In this paper, we start to system- 
atically establish a database of these and address the reliability 


1 Tubingen Model-Atmosphere WWW Interface, http : //astro . 

uni - tuebingen . de/~TMAW 
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Table 1 . Parameters of the HST DA standard stars (Gianninas et al. 
2011 ). 


Name 

WD no.“ 

T e ff 

log 9 



[K] 

[cm/s 2 * ] 

G191-B2B 

0501 +527 

60920 + 993 

7.55 ± 0.05 

GD71 

0549 + 158 

33 590 ± 483 

7.93 ± 0.05 

GD 153 

1254 + 223 

40320 + 626 

7.93 ± 0.05 

HZ 43 A 6 

1314 + 293 

56 800+ 1249‘- 

7.89 ± 0.07 


Notes. (a) WD numbers are from McCook & Sion (1999). (i,) HZ43A 
is only used in the UV because of contamination at longer wavelengths 
from its M-dwarf companion. (c) Beuermann et al. (2006, 2008) deter- 
mined T e g = 51 111 ± 660 K and log g- 7.90 ± 0.080. 


of state-of-the-art model-atmosphere spectra and the achievable 
limits in future flux calibration. 

White dwarfs (WDs) are ideal objects for the calibration 
of astronomical observations (Rauch 2012). They are relatively 
simple objects and their radiation is determined by fundamen- 
tal physics, e.g. their radius is defined by electron degeneracy. 
Moreover, they are nearby and their distance can be measured 
precisely, at least by the upcoming GAIA 2 mission (cf. Pancino 
et al. 2012, for a description of the GAIA spectrophotomet- 
ric standard stars survey). Most of the hot, hydrogen-rich WDs 
(spectral type DA) with T t s < 40000 K have virtually pure 
hydrogen atmospheres (gravitational settling), while the hotter 
WDs exhibit lines of heavier elements due to radiative levitation. 
WD spectral modeling requires adequate observations (WDs are 
intrinsically faint) and state-of-the-art theoretical model atmo- 
spheres that account for reliable physics and deviations from 
local thermodynamic equilibrium (LTE). 

The hot DA-type WD G191-B2B (BD+52°913) is, together 
with GD71, and GD 153, one of the primary flux reference 
standards for all absolute calibrations from 1000 to 25 000 A 
(Bohlin 2007). Recent results for their T e ff and loggr are sum- 
marized in Table 1. G191-B2B, the hottest and visually bright- 
est (my = 11.7228, van Leeuwen 2007) isolated WD (with a 
well known distance of 57.96 pc, Anderson & Francis 2012) 
of the sample, is ideal for panchromatic calibration from the ul- 
traviolet (UV) to the infrared (IR) wavelength range. However, 
due to its high T e ff and relatively low log g, radiative levitation 
competes against gravitational settling and holds trace elements 
in the photosphere and exhibits many weak metal lines (e.g. 
Barstow et al. 2003) in its observed UV spectrum. 

A variety of previous spectral analyses of G191-B2B 
(Table 2) had shown that it is difficult to determine its T t g- pre- 
cisely. Barstow et al. (1998) found that the metal content in 
the photosphere has a strong impact on the determined T e g. 
Teff = 60 920 K was found by the most recent analysis (Gianninas 
et al. 2011) who considered only C, N, and O (at solar abun- 
dances) in their models. The neglect of other metals calls for 
improved models with better metal opacities. The same may be 
true for HZ 43 A even if metals are below the detection limit of 
the available spectra (Table 1). 

Many abundance analyses were performed, most of them 
(e.g. Barstow et al. 2005), were based on previous T t g deter- 
minations from Balmer-line fits (cf. Table 2) and not from self- 
consistent fits to models with varying metal abundances. Lanz 
et al. (1996) measured He, C, N, O, Si, Fe, and Ni abundances, 
Holberg et al. (2003) determined abundances of C, N, O, Al, Si, 

2 http://www.esa.int/Our_Activities/Space_Science/ 

Gaia_overview 


Fe, and Ni and gave upper limits for Mg, Cr, Mn, and Co. The 
compiled abundances are listed in in Table 3. 

Based on a grid of state-of-the-art line-blanketed NLTE 
model atmospheres that include opacities of all identified metals, 
we perform a detailed spectral analysis. We describe the avail- 
able observations in Sect. 2, followed by a brief introduction to 
our model atmosphere code and the atomic data (Sect. 3). The 
spectral analysis is summarized in Sect. 4 and we end with our 
conclusions (Sect. 7). 

2. Observations 

2.1. FUSE data 

G191-B2B was observed many times over the course of the 
FUSE mission in the wavelength range 910 A-l 190 A, both for 
calibration purposes and for studies of the interstellar medium. 
For the present study, only observations obtained in the first 
eight months of the mission through the LWRS spectrograph 
aperture were analyzed. This time period included the major- 
ity of the LWRS exposure time obtained during the mission, and 
had the secondary benefit that the detectors had not yet suffered 
much degradation from gain sag. The observation IDs of the 
datasets were: M1010201, M1010202, M1030602, P1041203, 
and S3070101. 

Apart from a few quirks affecting the M 1010201 and 
Ml 030602 observations, which were among the first obtained 
during the mission, the quality of the data is excellent. No 
SiC data were obtained in observations M1010201 orM1030602 
as a result of channel mis-alignment. Otherwise: exposure-to- 
exposure variations in flux were typically well under 1%, in- 
dicating good channel alignment. The detector region used to 
record spectral image data for LiF2b was offset from the actual 
spectrum position in the M1010201 observation, so those spec- 
tra were discarded. The net exposure times were 33.3 ks for the 
SiC channels, 36 ks for LiF2b, and 40 ks for LiFl and LiF2a. 

No special processing was applied to data from individ- 
ual exposures. Raw data were processed with CalFUSE v3.2.3. 
Zero-point offsets in the wavelength scale were adjusted for each 
exposure by shifting each spectrum to coalign narrow interstel- 
lar absorption features. In order to assess the influence of geo- 
coronal airglow emission, spectra obtained during orbital day 
and night were combined separately. All the observations were 
obtained in spectral image (“histogram”) mode, so no informa- 
tion on photon arrival time was available within an exposure. 
However, the timeline table in the intermediate data files was 
examined for each exposure to determine the time spent in the 
“Day” and "Night” portions of the orbit. If the “Day” portion of 
such an exposure exceeded 15% of the total exposure duration, 
it was included with the other Day spectra. Because histogram 
mode exposures were short, most exposures were entirely Day 
or entirely Night. 

The individual exposures from all five observations were 
then combined to form composite Day and Night spectra for 
each channel. The Day and Night spectra were then compared 
at the locations of all the known airglow emission lines. If 
the Day spectra showed any excess flux in comparison to the 
Night spectra at those locations, the corresponding pixels in 
the Day spectra were flagged as bad and were not included in 
subsequent processing. Significant airglow emission during or- 
bital Day was seen for most observations at Hi Ly /? through 
Lyd, and Oi 44 988, 1027, 1028, 1039 A. Significant airglow 
was present during orbital night only at Ly/3; this affects the 
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interstellar-absorption profile but has no effect on our analysis 
of the photospheric spectrum. 

The final step was to combine the spectra from the four in- 
strument channels into a single composite spectrum. Because of 
residual distortions in the wavelength scale in each channel, ad- 
ditional shifts of localized regions of each spectrum were re- 
quired to coalign the spectra; such shifts were typically only 
one or two pixels. Bad pixels resulting from detector defects 
were flagged at this point and excluded from further process- 
ing. Finally, the spectra were resampled onto a common wave- 
length scale and combined, weighting by signal to noise on a 
pixel-by-pixel basis. 

The signal to noise of the final combined spectrum is limited 
by fixed-pattern noise in the detectors. The final spectrum has a 
minimum of roughly 20 000 counts per 0.0 1 3 A pixel in the con- 
tinuum, near the Lyman edge, and 60000-130000 counts per 
pixel long- ward of 1000 A. The effects of fixed-pattern noise 
are minimized by the fact that the positions of the spectra on 
the detectors varied during each observation, and by the fact 
that nearly every wavelength bin was sampled by at least two 
different detectors. 

2.2. HST data 

As described in detail by Bohlin & Gordon (in prep.), the 
HST/STIS low-dispersion flux calibration is derived from an 
ensemble match to the NLTE TLUSTY (version 203) model 
atmosphere SEDs for pure hydrogen (Hubeny & Lanz 1995). 
The models are for G191-B2B, GD 71. and GD 153. Originally, 
HZ 43 A was also used as a standard star but fell off the list of pri- 
mary flux standards because of an M star companion that con- 
taminates the STIS observations in the visible and IR (Bohlin 
et al. 2001). 

For the STIS echelle modes, the flux calibration is based only 
on the TLUSTY model for G191-B2B. The echelle absolute 
fluxes are less precise than for low dispersion because of the sin- 
gle model for the reference fluxes, because of imprecision in the 
matching of the separate echelle orders, and because the plethora 
of weak lines at the shorter wavelengths are missing in the refer- 
ence SED. However, the STIS echelle narrow metal line profiles 
are unaffected by uncertainties in the absolute fluxes. 

For the highest STIS resolution of ~3 kms -1 , there are 
two modes, namely E140H and E230H, which require sev- 
eral central wavelength settings for complete wavelength cov- 
erage from 1145-3145 A. Because G191-B2B is the primary 
STIS echelle calibration star with repeated observations, 105 ob- 
servations in the 072 x 072 aperture are available from the 
Mikulski Archive for Space Telescopes (MAST) 3 . Each spec- 
trum is resampled to a wavelength grid with a sampling inter- 
val corresponding to a resolving power of R = 2.3 x 10 5 and 
co-added. The number of individual observations at each wave- 
length point ranges from 4-44, while the total exposure time 
ranges from 6400-64 000 s at each point. The total counts in 
electrons at each point in the continuum are typically well above 
1000 and range up to over 10 000 from 1225-1400 A, where the 
statistical uncertainty is sometimes better than 1%. The high- 
dispersion echelle spectrum is available from the CALSPEC 4 
database along with the STIS low-dispersion data. 


3 http://archive.stsci.edu/ 

4 http: //www. stsci . edu/hst/observatory/cdbs/calspec . 

html 



AX/ A 

Fig. 2. Comparison of STIS and FUSE observations around Hi La-6 
with our final model. Thick (red in the online version) photospheric + 
ISM line absorption model with Nm = 2.2 x 10 18 cm -2 ; thin (blue) 
pure photospheric model; dashed (green. La only) IVm = 1.2 x 
10 18 cm -2 , 3.2 x 10 18 cm -2 . The locations of the D I blends are marked. 
L/3-6 are shifted in flux (0.7, 1.6, 2.7 x 10 -11 erg cm -2 s -1 A -1 ) for clar- 
ity. A reddening of E B -v = 0.0005 is applied following the law of 
Fitzpatrick (1999, withRy = 3.1). 

The photospheric radial velocity v r ^ = 22.1 ± 0.6 kms -1 
measured by Holberg et al. (1994) matches well our STIS 
observation. We adopt this value for our analysis. 

2.3. Interstellar line absorption and reddening 

The interstellar neutral hydrogen density Nm was determined 
from the comparison of our final model with the STIS and 
FUSE observations (Fig. 2). In all plots shown in this paper, we 
modeled the interstellar medium (ISM) line absorption (using 
Voigt line profiles) with WRPLOT 5 . The best match is found 
for log(AWcm 2 ) = 18.34+j}^. The Dl blends to Hi La-6 
are clearly visible and best reproduced at log(V Dl /cm 2 ) = 
13.54 4 ygg, i.e. D/H = 1.59 4 ggg x 10 -5 , Our values are in good 
agreement with those determined by Lemoine et al. (2002), 
log(AWcnr) = 18.18 ± 0.18 and D/H = 1.66 4 °® x 10 -5 (both 
with 2cr errors). 

Besides H I and D I, we identified interstellar lines of C II- IV, 
N I- II, Oi, Ain, Si II- ill, Pi- II, Si-n, and Fell in the FUSE 
and STIS spectra (Table 5). To identify pure photospheric lines 
that are contaminated by ISM lines, we modeled all of these and 
found that we need two distinct clouds with y ra d = 9 ± 1 km s -1 
and y la( j = 19 ± 1 km s -1 . This is well in agreement with the mea- 
surements of Sahu et al. (8.6 ± 1.7 kms -1 and 19.3 ±2.5 kms -1 
1999), who assigned the latter value to the local interstellar 
cloud. Dickinson et al. (2012b) measured 8.5 ± 0.18 kms -1 
and 19.3 ± 0.03 kms -1 . They unambiguously detected that the 
first cloud is of circumstellar origin. An additional third cloud 
with intermediate velocity like assumed by Vidal-Madjar et al. 
(1998, u ra d = 8.2, 13.2, 20.3 (±0.8 kms -1 ) is not necessary for 
our modeling (Fig. 3, top). 

Interestingly, we find additional weak absorptions of 
OiT 1302.163 A and Ni A 1 199.550, 1200.223, 1200.710 A 
at ty ac j of -26.3 kms -1 and -26.1, kms -1 , respectively. These 

5 http : //www. astro .physik. uni -potsdam. de/~htodt/ 
wrplot/index . html 
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Fig. 3. STIS spectrum around the interstellar absorption lines 
Ol A 1302.163 A (top) and Nl AA 1 199.550, 1200.223, 1200.710 A 
(bottom) compared with the synthetic spectrum of our final model 
where the ISM lines are included. The labels give the radial velocities 
(in km s -1 ) that are applied. 



log (X / A) 

Fig. 4. Comparison of FUSE and HST (STIS and NICMOS) observa- 
tions with our final model. A reddening of E B _ V = 0.0005 is applied. 
The low -resolution (LR) STIS+NICMOS observation vanishes behind 
the model SED due to the line width. Therefore, we plotted the ob- 
served spectrum twice, one shifted by A log J\ = -0.2 for clarity. U, R , I 
(Landolt & Uomoto 2007), B, V (H0g et al. 2000), J, H, and K (Cutri 
et al. 2003) fluxes (converted from brightnesses using values given by 
Heber et al. 2002) are shown for comparison. 


velocities are reminiscent of the expansion velocity of a plane- 
tary nebula shell (e.g.. Kwok et al. 1978), that for a stellar mass 
of M = 0.555 M q (Sect. 4.6) must have been ejected more than 
500 000 years ago (Renedo et al. 2010). Its recombined, neutral 
gas, however, is still in the line of sight. 

From the low interstellar Nu i density, we expect a low inter- 
stellar reddening. The Galactic reddening law of Groenewegen 
& Lamers (1989), log (Nm/E B -v) = 21.58 ± 0.10, predicts 
0.0003 < Eb-v Si 0.0007. Figure 4 shows a comparison of 
observations and synthetic spectrum from the far UV (FUV) to 
the IR. We find E B . V = 0.0005 ± 0.0005. 


3. Model atmospheres and atomic data 

Table 2 demonstrates clearly that a panchromatic analysis from 
the EUV to the optical is inevitable for accurate results on pho- 
tospheric parameters. Moreover, NLTE modeling is mandatory 
to calculate a reliable synthetic spectrum. 

Lanz et al. (1996) presented the first NLTE model (Table 2) 
that reproduced the observed spectrum from the EUV to the op- 
tical wavelength range. Barstow & Hubeny (1998) introduced 
then a stratified H+He envelope including heavier metals in 
their models to improve the match to the observed flux below 
the Hell absorption threshold (T 5 228 A). Later analyses 
had shown that there is further evidence for a stratification in 
G191-B2B’s photosphere. Vennes et al. (2000) closely exam- 
ined Feige24 that, compared with G191-B2B, has similar at- 
mospheric parameters and an almost identical abundance pat- 
tern. They found that the O IV / O V ionization equilibrium is 
overcorrected by -0.8 dex in their NLTE model. They concluded 
that this might reveal an inhomogeneous vertical stratification of 
oxygen in both stars. A later analysis of both stars (Vennes & 
Lanz 2001) showed that the average heavy-metal abundance in 
Feige24 is 0.17 dex larger compared to the cooler and, hence, 
older G191-B2B (same loggr). Thus, the abundance pattern is 
determined by the same processes in both stars and the authors 
assumed that selective radiative pressure and gravity are in diffu- 
sive equilibrium. This was proven by Dreizler & Wolff (1999). 
They used self-consistent diffusion models (Table 2) that were 
able to reproduce the observed flux for A £ 228 A without ad- 
ditional absorbers or mechanisms. However, problems remained 
with the fit to the UV lines. 

Now, our strategy to proceed with the analysis is threefold. 
We start with chemically homogeneous models to find the model 
that reproduces best the continuum slope and the spectral lines 
from the FUV to the optical (Sect. 4). In an intermediate step, 
we will then apply the depth-dependent abundance profiles cal- 
culated by Dreizler & Wolff (1999) to our final homogeneous 
model to investigate the impact of chemical stratification on 
the emergent spectrum (Sect. 4.3). In the last step, a diffusion 
model is calculated and compared with the homogeneous model 
(Sect. 4.4). 


4. Spectral analysis and results 

The metal-line blanketed NLTE model atmospheres for our anal- 
ysis were calculated with the state-of-the-art Tubingen NLTE 
model-atmosphere package (TMAP 6 , Werner et al. 2003), which 
can consider opacities of all elements from H to Ni and be- 
yond (Rauch 1997, 2003; Werner et al. 2012; Rauch et al. 2012). 
TMAP was successfully used for the spectral analysis of hot, 
compact stars (e.g. Rauch et al. 2007; Wassermann et al. 2010; 
Ziegler et al. 2012). 

Our models assume plane-parallel geometry and are in hy- 
drostatic and radiative equilibrium. Opacities of all species for 
which spectral lines are identified, namely H, He, C, N, O, Al, 
Si, P, S, Ca, Sc, Ti, V, Cr, Mn Fe, Co, Ni, Zn, Ge, and Sn, 
were considered in the model-atmosphere calculations. For all 
elements, we account for level dissolution (pressure ionization) 
following Hummer & Mihalas (1988) and Hubeny et al. (1994). 
Figure 5 demonstrates that our H I model ion (Table 4) includes 
all levels that are relevant in the line-forming region -4.5 < 
log [w/(g/cm 2 )]. All model atmospheres presented here cover 


6 http : //astro . uni - tuebingen . de/~TMAP 
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Fig. 5. Occupation probabilities of the H I levels with principal quantum 
numbers n = 1-14 in our final model. 


column densities m of -7.6 < logm < 3.2 (cf. Beuermann et al. 
2006) represented by 90 depth points. 

The model-atoms and respective absorption cross-sections 
for Ca-Ni were calculated via the recently registered VO service 
TIRO that uses Kurucz’s atomic data 7 8 9 10 and line lists (Kurucz 
1991,2009, 2011). 

Ca, Sc, Ti, V, Cr, Mn, and Co lines are not identified. These 
were merged into a generic model atom (Rauch & Deetjen 2003) 
with fixed solar abundance ratios. Then, we performed test cal- 
culations and adjusted the abundance to a value (1.78 x 10 -6 by 
mass, the solar value is 9.93 x 10 -5 ) where all of its lines just 
fade in the noise of the observed spectra. All other model atoms 
were constructed from data retrieved from the public Tubingen 
model-atom database TM AD 1 ' . 

In total, we considered 1038 atomic levels in NLTE com- 
bined with 4646 line transitions (for the number of individual 
iron-group lines, see Table 4) in the model-atmosphere calcula- 
tions with 53 203 frequency points within 1 x 10 12 Hz < v < 
3 x 10 17 Hz. For the emergent spectra (100 A < A < 400 000 A, 
686 196 wavelength points), we account for fine-structure split- 
ting and used 1585 NLTE levels and 9721 respective line tran- 
sitions. The model-atom statistics are summarized in Table 4. 
Figure 6 shows the ionization fractions of all elements in our fi- 
nal model. It may be interesting to note that a single model atmo- 
sphere needs about one week to converge, i.e. the absolute values 
of all relative corrections are below 10 -4 , on a 64 bit, 2.66 GHz 
compute core with 8 GB memory. 

For the calculation of synthetic Hi line profiles, we use 
Stark line-broadening tables provided by Tremblay & Bergeron 
(2009). For those lines, where no broadening tables are available, 
TMAP uses an approximate formula, as described in Ziegler 
et al. (2012, Eqs. ( 1 )— (5)). 

We started with a model with T e ff = 60 920 K and log g = 7.55 
(the values of Gianninas et al. 2011) and the element abun- 
dances from Table 3. Next, we adjusted these abundances to best 
reproduce the respective spectral lines. We then calculated an 
extended grid of 234 model atmospheres (48 000 K < T e g < 
68 000 K in steps of <1000K and 7.35 < log g < 7.75 in 
steps of 0.05 (some of the hotter models are calculated only 
for log g < 7.60). For this grid, we extensively used compute 
resources of the bwGRiD 1 " in addition. Although this highly 
speeded up the model-grid calculation, the wide parameter range 

7 http : //astro . uni - tuebingen . de/~TIR0 

8 http : //kurucz . harvard . edu/atoms . html 

9 http : //astro . uni - tuebingen . de/~TMAD 

10 http : //www . bw- grid . de/en/ the-bwgrid/ 


and the large number (15) of parameters to adjust simultaneously 
did not allow us to take a statistical approach in the spectral anal- 
ysis (x 2 method like e.g. in Gianninas et al. 201 1) on a reason- 
able time scale. We therefore need to rely upon our “^-by-eye” 
methods. All SEDs that were calculated for this analysis are 
available via TheoSSA 11 . 

In a first analysis step, we will determine log g based on UV 
and optical observations. Then, we will determine T e g precisely 
based on ionization equilibria of the metals which are sensitive 
indicators. Subsequently, we will adjust the abundances again 
and verify our T e g and log g results. 

4.1. Surface gravity and effective temperature 

The dependency of the synthetic flux level on T e ff and log <7 
for fixed abundances is demonstrated in Fig. 7, where we com- 
pare the observed and synthetic fluxes in the FUV. In the top 
panel, it is obvious that at a constant T e u = 60 920 K, a log g 
higher by 0.2 dex than log <7 = 7.55 measured by Gianninas et al. 
(2011) is necessary to reproduce the Lyman-line decrement. For 
a fixed log <7 = 7.55 (middle panel), a lower Teff (A T e ff = 6000 K) 
improves the agreement between model and observation. The 
bottom panel shows that at values within the (statistical) error 
ranges from the Hi Balmer-line analysis, T e ff = 60000K and 
log <7 = 7. 60 (cf. Table 2 Gianninas et al. 2011), a good agree- 
ment for both, line profiles and decrement, is achieved. 

This was not expected from the outset although Barstow 
et al. (1998) found a relatively good agreement of T e ff and log <7 
from H I Lyman and Balmer lines in the heavy-metal rich mod- 
els (Table 2). The later analysis by Barstow et al. (2001, Table 2) 
shows strong deviations in log g between optical and FUV anal- 
yses. Figure 8 shows a comparison of synthetic H I Balmer line 
profiles with optical observations. The deviation between the 
T e ff = 60 920 K / log <7 = 7.55 and the 60 000/7.60 ones is minor. 

In addition to the low-resolution ( R ~ 500) optical spec- 
trum, medium-resolution ( R ~ 6000) observations of Ha and 
HyS are shown is Fig. 9. The agreement among the STIS low and 
medium resolution is excellent. While the model absorption is a 
bit weak as shown in the lower plot of Fig. 9, the central NLTE 
emission reversal agrees well with the observations. 

Although we cannot reproduce Ha and H )3 in detail in the 
medium-resolution spectrum, this has no significant influence 
on our determination of T e tf and log <7 because the higher mem- 
bers of the Hi Balmer series form much deeper in the atmo- 
sphere where the influence of metal opacities is less important 
(cf. Napiwotzki & Rauch 1994). We adopt log g = 7.60. 

In the next step of this analysis, we evaluate ioniza- 
tion equilibria of metals that exhibit lines of successive ion- 
ization stages. Figures 10-12 show some strategic lines for 
this T e ff determination. In total, we can use eight elements and 
lines of Cm- IV, Nm-v, Sim-iv, Piv-v, Siv-v, Feiv-vi, 
Ni IV- VI, and Ge IV- V. T e g = 60 000 ± 2000 K reproduces well 
all these equilibria simultaneously. For our further analysis, we 
adopt T e ff = 60 000 K. 

4.2. Photospheric abundances 

In the following, we use logarithmic mass fractions for all abun- 
dance values, if not otherwise mentioned. Previously determined 
abundances and respective references are summarized in Table 3. 
In the following, we will briefly mention the strategic lines for 


11 http://dc.g-vo.org/theossa 
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Fig. 7. Section of the FUSE observation compared with our model fluxes with different T e g and log g. In the top and middle panels , the synthetic 
fluxes are normalized to the observed flux at 1000 A and in the bottom panel to the observed K magnitude (see Fig. 4). E B -v and Am are applied 
using our results from Sect. 2.3. 


the abundance determinations and note abnormalities for an el- 
ement selection only. Most of the identified metals exhibit lines 
of at least two subsequent ionization stages and some of these 
lines were already used for the determination of T e ff (Sect. 4.1). 
The abundances were then adjusted to achieve best line fits. 
Two large plots (Figs. A.l, A. 2, German DIN A0 size) are pro- 
vided in the online material that show a comparison of our final 
model with the observation in the FUSE and STIS wavelength 
ranges (in total 911-1750 A). They include all line identifica- 
tions (FUSE/STIS wavelength range), e.g. 2/421 Feiv, 144/815 
Fev, 1/52 Fe VI. 1/236 Niiv, 13/690 Ni V, and 9/43 Ni VI lines. 
These numbers are much higher than those of Preval et al. (2013, 
106 Fev and 44 Ni V lines in the STIS wavelength range). The 
recent work of Berengut et al. (2013) to employ G191-B2B as 
a stellar laboratory to determine the fine-structure constant is 
based on the latter list and may, thus, not fully exploit capacity 
of all the available STIS spectra of G191-B2B. 

Our line identifications are also summarized in Table 5, 
whereas Table 6 gives a list of the strongest unidentified lines. 


4.2.1. Helium 

The first analyses revealed only upper limits for the He abun- 
dance, e.g. He < -3.1 and <4.1 (Vennes et al. 1996; Gunderson 
et al. 2001, respectively). Cruddace et al. (2002) determined 
He = -4.2 ±0.1 using high-resolution EUV spectroscopy. An 
attempt to identify and measure He II Lyman lines (n - n' = 1-4, 
1-5) with J-PEX 12 (Barstow et al. 2005) was not successful. Our 
models show that Hen A 1640 A (2-3) should be clearly visible 
at He = -3.7 and -4.2 and disappears in the noise of the observa- 
tion only at about He < -4.7 (Fig. 13). We adopt this upper-limit 
value for our models. 

4.2.2. Carbon, nitrogen, and oxygen 

C III and C IV lines are visible in the observation. C hi A 977.02 A 
and C iv AA 1548.20, 1550.77 A have strong ISM blends. In case 


12 Joint Astrophysical Plasmadynamic Experiment 
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Fig. 8. Synthetic Hi Balmer-line profiles compared with the STIS LR 
observation. The observation is shifted by factors of 0.921, 0.323, 
0.232, 0.202, 0.193, and 0.188 for Her - Hf, respectively, to fit into 
this plot. Red, full line: r e (f = 60000K and log g = 7.60; blue, dashed: 
T tS = 60 920 K and logg = 7.55. 



Ml A 


Fig. 9. STIS Her and Hj3 low-resolution (dashed, blue) and medium- 
resolution (gray lines) observations (labeled with medium-resolution 
grating/central wavelength in A). Because of uncertainty in the flux 
calibration, the medium-resolution data are normalized to the low- 
resolution flux. The flux around H/3 in the top plot is multiplied by a 
factor of 0.35. The red lines in the upper two plots are the medium- 
resolution ( R a 6000) spectra degraded to the low resolution (R » 500) 
and agree with the low resolution (blue dash) within the uncertainty of 
the R = 500 resolution. The lower two plots are shifted down by 0.035 
and 0.07 x 10 -13 flux units, respectively. In the lowest plot, the model is 
overplotted in red after smoothing to the medium resolution. While the 
model Her absorption is somewhat too weak, the central emission agrees 
with the observation within the uncertainty of the resolution (insert). 

of the latter, the photospheric component can be separated and 
modeled (Fig. 10). At C = -5.15, lines of both ions are well 
reproduced. 

N III- V lines are found in the observation, they are all well 
matched at N = -5.58 (Fig. 10). 

Vennes et al. (2000) encountered deviations between oxy- 
gen abundances determined from OlV and Ov lines in an 


Table 6. Wavelengths (in A) of unidentified strong {W,\ > 10 mA), 
likely photospheric lines in the FUSE and STIS observations. 


FUSE 

989.11 

1029.44 

1133.04 

1142.91 


STIS 

1157.60 

1173.29 

1183.82 

1201.51 

1283.52 

1158.08 

1174.34 

1186.08 

1201.81 

1306.04 

1158.76 

1176.11 

1186.27 

1202.43 

1331.19 

1165.31 

1176.52 

1186.59 

1204.47 

1385.32 

1166.80 

1176.64 

1187.70 

1227.55 

1389.89 

1171.12 

1176.87 

1190.15 

1253.67 

1398.21 

1171.19 

1176.98 

1192.01 

1253.87 

1411.46 

1171.45 

1177.06 

1194.16 

1258.81 

1516.68 

1172.18 

1178.68 

1198.15 

1260.48 

1520.64 

1172.34 

1182.00 

1201.30 

1272.98 

1525.32 

1173.22 

1183.37 

1201.47 

1281.37 

1538.94 



XI A 


Fig. 13. Synthetic spectrum around Hell A 1640.42 A compared with 
the STIS observation. He abundances are green, dashed: -10.0, red, 
thick: -4.7, blue, thin: -4.2, blue, dashed: -3.1. The insert shows the 
region Ad = ±0.4 A around the Hell line. For comparison, the obser- 
vation was smoothed with a low-pass filter (Savitzky & Golay 1964, 
n = 15, m = 4). 


analysis of Feige 24. The O V abundance was 0.5 dex higher 
in their LTE model approach. In their NLTE models, they 
found that the O IV / O V ionization equilibrium was overcor- 
rected by -0.8 dex. They suggested an inhomogeneous strati- 
fication of O in the atmosphere. Vennes & Lanz (2001) dis- 
covered that a similar problem exists in G191-B2B, with an 
overcorrection of -0.6 dex. Consequently they assumed that in 
both stars, the interplay between selective radiation pressure and 
gravity in diffusive equilibrium are the key processes for this 
phenomenon. Figure 14 shows the same deviation in our mod- 
els. While Otv AA 1338.615, 1342.990, 1343.526 A are well fit- 
ted at O = -4.72, Ovl 1371.296 A is apparently much stronger 
than observed. It is matched with an O abundance that is reduced 
by -0.4 dex. 

In the FUSE observation, only the short wavelength com- 
ponent of the Ovidd 1031.912, 1037.614 A resonance doublet 
is detectable. The unexpected weakness of this doublet was al- 
ready reported by Oegerle et al. (2005). Dickinson et al. (2012b) 
verified that it stems from the photosphere. The Ovi reso- 
nance doublet in our models is even stronger, compared to 
OlV and Ov lines, requiring a reduction of the O abundance 
by about -1.5 dex (Fig. 14). Dickinson et al. (2012a) encoun- 
tered a similar problem with enigmatically deep line profiles of 
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Fig. 14. Most prominent O IV, O V, and O VI lines in the STIS and FUSE 
observations compared with our theoretical line profiles calculated with 
three O abundances (O = -4.72, -5.12, -6.24, from top to bottom). 



Fig. 15. Theoretical A1 III line profiles calculated from our final model 
compared with the STIS observation. 


the N V resonance doublet in their models. We revisit the prob- 
lem with the oxygen abundances derived from different ioniza- 
tion stages in Sects. 4.3 and 4.4 in detail. 

4.2.3. Aluminum, silicon, phosphorus, and sulfur 

Holberg et al. (1998) identified the Aim AA 1854.72, 1862.79 A 
resonance doublet in the IUE NEWSIPS SWP Echelle Data 
Set 13 , and Holberg et al. (2003) measured Al = -5.08. We could 
newly identify some other Aim lines. We derive Al = -4.95, 
well in agreement with the Holberg et al. (2003) value (Fig. 15). 

Si III- IV, P IV- V, and S V- VI lines are identified. We deter- 
mine Si = -4.30, P = -5.81, and S = -5.24 (Fig. 10). 


13 http: //vega . lpl . arizona.edu/newsips/ 


4.2.4. Iron-group elements 

Many hundreds of lines of FelV-Vl and NilV-vi are iden- 
tified (Table 5). They are best reproduced at Fe = -3.30 
and Ni = -4.45. Note that the Ni/Fe abundance ratio is 
about 25% higher than the solar ratio. Some of these lines are 
shown in various figures in this paper, please have a look at 
the two large online figures that show the complete FUSE and 
STIS wavelength ranges. An animation of STIS wavelength 
range can been seen at http://astro.uni-tuebingen.de/ 
~rauch/A0_E14®H_SW. gif as well. 

In Fig. 15 three lines of CrlV and one of Co IV are visible in 
the synthetic spectrum of our final model. These are weak and 
comparable to the noise of the observation. Although one may be 
tempted to believe the presence of Criv A 1863.075 A, we take 
this as a hint that a log mass fraction of -5.75 for the generic 
model atom is reasonable and adopt this as an upper limit for our 
analysis. This is, within the error limits, in agreement with the 
upper limits for Cr, Mn, and Co of about -6.2 that were found 
by Holberg et al. (2003). 

Preval et al. (2013) suggested that the unidentified line 
at 1272.98 A is a V IV line. Since many other V IV lines 
with much stronger log gf values ( g is the statistical weight 
of the lower atomic level and / is the oscillator strength 
of the line transition) from Kurucz’s POS line lists (with 
good wavelengths) are not present in the spectrum, e.g. 
Viv/tT 1355.127, 1419.577, 1426.647 A (all more than ten 
times higher log gf) therefore this identification appears to be 
very unlikely. 

We mention here that we find deviations between Kurucz’s 
POS wavelengths and the observation of up to 0.05 A. In addi- 
tion, Fig. 13 shows that the strengths of Feiv A 1640.042 A and 
Feiv A 1640.155 A in the model are the opposite way around in 
the observation. 

Figure 16 shows a comparison of models (calculated with 
Kurucz’s POS lines) in the FUSE and STIS wavelength ranges 
where in each case the abundance of an individual element X in 
the construction of the generic (Ca, Sc, Ti, V, Cr, Mn, Co) model 
atom is increased by a factor of ten. Values higher that 1 in the 
flux ratio indicate stronger lines of element X. 

E.g. the case of Ti, two lines are much stronger than all oth- 
ers, Tiiv AA 1451.739, 1467.343 A. They are not identified in the 
observation but at the ten times increased abundance they are 
clearly visible in the model. The same is valid for Cr, where 
Cr iv AA 1332.415 A and Crvi AA 1417.660 A are the strongest 
lines in our models (Fig. 16), and for Mn and Co as well. This 
allows us to establish upper abundance limits of about 10% solar 
for Ti, Cr, Mn, and Co (cf. the beginning of Sect. 4). 

4.2.5. Zinc, germanium, and tin 

21 Zniv lines are newly identified in the STIS observation. 
These are almost all that are listed in the NIST 14 database with 
relative intensities higher than 100. Since no individual calcula- 
tions for Zn IV transition probabilities are available, we adapted 
those of the isoelectronic Ge VI (Rauch et al. 2012). In Fig. 17, 
we show nine of them with NIST relative intensities of 200. All 
their theoretical line profiles are reproduced at Zn = -4.89. 

For Ge, we used the same model atom as Rauch et al. (2012) 
and determined Ge = -5.49 (Fig. 10). 


14 http : //physics .ni st . gov/PhysRefData/ASD/lines_form. 
html 
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Fig. 16. Flux ratio of our final model and a model with ten times in- 
creased abundance of element X (Ca to Co, from top to bottom). The 
location of V IV 41272.98 A is marked. 

We constructed a relatively small Sn model atom. The only 
lines for which reliable oscillator strengths are available are the 
Snm and Sniv resonance lines (Morton 2000). For all other al- 
lowed transitions, we follow Werner et al. (2012) and set / = 1 . 
We used the Sniv A 1314.537 A resonance line, like Vennes et al. 
(2005), to measure the abundance of Sn = -6.45. 

4.2.6. Summary of results with chemically homogeneous 
models 

We can reproduce the entire ultraviolet spectrum of G191-B2B 
with our chemically homogeneous NLTE models, with the 
exception of the Oiv/vi lines which are obviously af- 
fected by O stratification effects. Current diffusion models 
yield poor fits to the metal lines (Dreizler & Wolff 1999). 
T e ff = 60 000 ± 2000 K and log ry = 7.60 ± 0.05 were determined 
within small error limits. They are in agreement with Gianninas 
et al. (2011, r e ff = 60 920 ± 993 K, log g = 7.55 ± 0.05). We do 
not encounter problems in modeling Hi Lyman and Balmer 
lines simultaneously with the same T e ff and log g like found by 
Barstow et al. (2001, see Table 2). 


We can determine all abundances with error limits of 0.2 dex. 
In case of Zn, where we adopt Ge VI /-values, we estimate that 
the error is 0.3 dex. Our C, N, O, Al, Si, Fe, and Ni abundances 
(Fig. 18) agree, within error limits, with those of Vennes et al. 
(1996); Holberg et al. (2003); Vennes et al. (2005). Our val- 
ues are in general slightly higher. One reason may be the about 
6000 K higher T e g of our final model. The stellar parameters are 
summarized in Table 8. 

The abundances of all elements but Fe predicted by Chayer 
et al. (1995) for a DA-type WD differ strongly from those that 
we determined (Fig. 18). 

4.3. Test of the diffusion impact 

In a first step, we simply applied the abundance profiles provided 
by Dreizler & Wolff (1999) to the occupation numbers of He, C, 
N, O, Si, and Ni in our final model. Figure 19 (top panel) shows 
that this gives a good agreement with O V while O IV is now too 
weak. The O VI lines appears even stronger, strengthening the 
discrepancy. Since the atmospheric structure was kept fixed in 
this test, we expected that, if at all, only a self-consistent dif- 
fusion model is able to reproduce the observed OlV-Vl lines 
simultaneously. 

4.4. A self-consistent diffusion model 

We used the NGRT 1 ' code (Dreizler & Wolff 1999; Schuh et al. 
2002) to calculate diffusion models with the same element com- 
position and model atoms like our homogeneous TMAP mod- 
els. The first model shows a strongly increased abundance of 
the generic model atom that combines Ca, Sc, Ti, V, Cr, Mn, 
and Co (Sect. 4) and, hence, much too strong lines of the con- 
sidered elements. The reason is that the IrOnlc code (Rauch & 
Deetjen 2003) calculates a mean atomic weight for the generic 
atom following 


where r, is the relative mass-fraction (with respect to r\ = 1) 
and A; the atomic weight of element i. The artificially increased 
number of lines of a single generic element strongly increases its 
radiative levitation. Flux blocking by the generic element then 
leads to stronger gravitational settling of other elements, e.g. 
Sn had an abundance below 10 -17 throughout the model atmo- 
sphere. The other elements showed abundances that were partly 
more than one dex below those of our homogeneous model. 
Since we did not want to neglect all opacities of the generic 
atom, we changed its atomic weight to 

n 

tIig = 'y, n ■ Ai. (2) 

1=1 

Now, the stratified NGRT models yields depth dependent abun- 
dances (Fig. 20) that are closer to those of our homogeneous 
model, especially Sn appears at a realistic value. In case of He, 
C, N, O, Si, and Ni the abundance profiles are similar to those 
of (Dreizler & Wolff 1999). The changed atmospheric struc- 
ture is shown in Fig. 21. It is interesting to note that most of 
the lines and all continua are formed at logm > -3 (Fig. 22) 
while deviations in the temperature structure are noticeable only 
outside of this region. The resulting spectrum (Fig. 23) of the 

15 New generation radiative transport. 
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Fig. 17. Theoretical line profiles of the strongest ZnlV lines in the STIS wavelength range compared with the observation. 



Fig. 18. Top : photospheric abundances of G191-B2B (red stars) com- 
pared with solar values (Asplund et al. 2009). [X] denotes log (mass 
fraction/solar mass fraction) of element X. The dashed, green line 
shows the solar ratio. The arrows indicate upper limits. The cyan dia- 
monds (Holberg et al. 2003), triangles (Vennes et al. 1996), and tridents 
(Vennes et al. 2005) are previously determined values. Bottom : com- 
parison of our abundance number ratios (red stars) with predictions of 
diffusion calculations for DA-type (blue squares) WDs (Chayer et al. 
1995) with T e ff = 60000 K and log g = 7.5. 


stratified model is, compared with the homogeneous model, no 
improvement. While Fe V lines match the observation at about 
T e ff = 55 000 K, it can be extrapolated that Fe IV lines are much 
too strong for T e g < 70 000 K. Ni IV and Ni V lines are much 
too strong because the Ni abundance is enhanced (Fig. 20) in 
the line-forming regions and can, thus, not be used for a T e g 
estimate. 

In the stratified models, the OlV and Ov lines are now 
mucher stronger than observed and Ovi appears at the same 
strength that resulted from our diffusion test (Sect. 4.3). 
The O abundance profile (Fig. 20) shows a strong increase 
for log m < - 4. Only by the introduction of an artificial abun- 
dance reduction by a factor of m/1585 for logm < -3.2, we 
achieve an acceptable agreement of O V and O VI (Fig. 19). O IV 
is still slightly too strong because the abundance and, thus, the 
lines (including a blend at O IV) of the generic iron-group atom 



AX/ A 


Fig. 19. Same as Fig. 14. Top: dashed, green TMAP (=chemically 
homogeneous) model with O = -4.72; thin, blue TMAP model with 
abundance profiles from Dreizler & Wolff (1999), thick, red NGRT 
(=diffusion) model. Bottom: thick, red NGRT model with an artificially 
reduced O abundance in the outer atmosphere. Note that in the NGRT 
models, the line strengths of the generic iron-group element (see text) 
are overestimated. 


(Sect. 3) are overestimated by the NGRT model (Fig. 20). Based 
on this numerical exercise, it may be speculated that a weak stel- 
lar wind or an other, unknown process that is not considered by 
NGRT is responsible for the lower oxygen abundances in the 
outer atmosphere. 

We can conclude two things. A generic model atom is obvi- 
ously not suited for a diffusion calculation due to the strongly 
enhanced number of lines for a single atom in the modeling 
process. The NGRT diffusion models yield partly too low abun- 
dances in the line-forming regions and, thus, cannot reproduce 
the metal line properly. An additional, weak wind may be nec- 
essary to increase the metal abundances in the line-forming 
regions. 

4.5. The extreme-ultraviolet spectrum 

The inability to model the EUVE spectrum with chemically ho- 
mogeneous atmospheres (Holberg et al. 1989) was the reason 
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Fig. 20. Abundance profiles in our diffusion model (Fuff = 60 000 K, 
log g = 7.60). Short-dashed (blue) lines: unrealistically high abundance 
of the generic iron-group element (IG, see text), thick (red) lines: re- 
duced IG abundance, horizontal long-dashed, thin (cyan) lines: the 
abundances in our final homogeneous model. In the O panel, the thick 
dashed (green) line shows our modified O-abundance profile (see text). 


to investigate stratified photospheres (e.g. Koester 1991). Lanz 
et al. (1996) demonstrated, that it is possible to consistently 
match the optical, UV, and EUV data with homogeneous NLTE 
models with the same T t s and chemical composition. 

We calculated EUV spectra from our model grid with 
193 584 frequency points within 100 A < A < 930 A, and 
Kurucz’s LIN line lists (theoretical and laboratory measured 
lines, in total 8 135405 lines of Ca-Ni in our wavelength 
interval, Kurucz 2009). These spectra were processed with the 
recently registered VO tool TEUV lf ' that corrects synthetic stel- 
lar fluxes for interstellar absorption below 911 A. It simulates 
radiative bound-free absorption of the lowest ionization states of 
H, He, C, N, and O using Opacity Project data (Seaton et al. 
1994). Two interstellar components with different radial and 
turbulent velocities, temperatures, and column densities can be 


16 http : //astro . uni- tuebingen . de/~TEUV 



Fig. 21. Temperature structures (T eff = 60 000 K, log g = 7.60) of a pure 
H model (thin, black), our final homogeneous model (thick, red), and 
a diffusion model (dashed, blue). The formation depths of the line 
cores of the lowest members of the Hi Lyman and Balmer series, the 
C IV, N V, O VI, Si IV, and S VI resonance doublets, and our strategic 
Oiv AA 1338.634, 1343.022, 1343.526 A and Ov 4 1371.296 A lines 
are marked. 



considered. Figure 24 shows the comparison of synthetic and 
observed EUV spectra. Our synthetic spectra were normalized 
to the measured FUSE flux of 1.347 x 10 _11 ergcm _2 s _1 A -1 at 
920 A. Then, the interstellar column densities are adjusted, to 
match the EUVE flux Nm for 530 A, (Vnei for 470 A, and Nhcu 
for 220 A. Since our models do not reproduce the measured flux 
between 250 A and the He II ground state threshold, V He n is not 
reliable. Table 7 shows the applied N m and V He i values com- 
pared with the literature values. Our N m values, necessary to 
match the EUVE flux level, are about a factor of two higher 
than log (Vhi / cm 2 ) = 18.34^°p that we determined previ- 
ously from Hi Lyman-line fits (Sect. 2.3). logiVNi = 13.87 
and log Ncn = 14.86 were adopted from Lemoine et al. (2002). 

The overall agreement of our homogeneous models with 
T e ff = 60 000 K at wavelengths A > 250 A is very good, espe- 
cially the interval 360 A < A < 450 A is excellently matched 
in detail. Models with T e g = 65 000 K and T e g = 55 000 K yield 
much too high and too low fluxes, respectively. At A < 250 A the 
theoretical flux is too high in all models, even at r e ff = 55 000 K. 
A stratified model (Fig. 24) with Teg- = 60 000 K fails to repro- 
duce the flux between 250 A < A < 420 A and has a too-high 
flux at A < 200 A. 

Both, our homogeneous and our stratified models, fail to re- 
produce the entire EUV spectrum of G191-B2B. It seems likely 
that there may be some stratification in the atmosphere but we 
don’t yet know how to distribute the various atomic species with 
depth. This is a challenge for theorists. 
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Fig. 24. Top : comparison of three synthetic spectra (r eff = 55 000, 
60 000, and 65 000 K) of chemically homogeneous models with the 
EUVE observation. The wavelengths of ground-state absorption thresh- 
olds of He ions are indicated. Bottom', comparison of a stratified model 
T e ff = 60 000 K with the observation. 


Table 7. Logarithmic ISM column densities for homogeneous models 
with different T eS to match the EUVE flux level of G191-B2B. 


Teff/X 

N m 

A^Hel 


Homogeneous (TMAP) 


55 000 

18.53 

17.45 


60000 

18.59 

17.45 


65 000 

18.64 

17.45 


Stratified (NGRT) 



60000 

18.60 

17.45 


Literature 




59250 

18.23 

17.16 

Kimble et al. (1993, HUT) 

54000 

18.27 

17.16 

Dupuis et al. (1995, EUVE) 

55 200 

18.32 

17.26 

Lanz et al. (1996, EUVE) 

56000 

18.32 

17.15 

Dreizler & Wolff (1999, EUVE) 

53 000 

18.28 

17.16 

Vennes & Lanz (2001, EUVE) 

54000 

18.33 

17.34 

Cruddace et al. (2002, J-PEX) 


4.6. Mass and distance 

A stellar mass of M = 0.555^'^ M 0 and a luminosity of 
\og(L/Lo) = 0.63^34 are determined by comparison with evo- 
lutionary models (Fig. 25) for old white dwarfs (metallicity 
z = 0 . 001 ). 

We calculated the spectroscopic distance following the flux 
calibration of Heber et al. (1984b) for T e ff = 5454 A, 

d[ pc] = 7.1 1 x 10" 4 ■ ^H v ■ M x io o - 4m vo- lo g0, (3) 

with my, = my - 2.175c, c = \A1Eb-v, and the Eddington 
flux H v = 1.109 x 10 -3 ergcm _2 s _1 Hz -1 at W of our final 
model atmosphere. We used Eb-v = 0.0005 ± 0.0005 (Sect. 2.3), 
M = 0.555+“ “ 3 M 0 , and m v = 1 1.7228 ± 0.0082 (van Leeuwen 
2007), and derived a distance of d = 62 ± 4 pc and a height above 



Fig. 25. Location of G191— B2B in the log /..fj-log g plane (the ellipse 
indicates the errors of our analysis) compared with evolutionary tracks 
for hydrogen-rich white dwarfs (Renedo et al. 2010) labeled with the 
respective stellar masses (in M 0 ). 

Table 8. Parameters of G191-B2B as derived by our analysis. 


Teff/ K 
log b/cm/s 2 ) 

60 000 ±2000 
7.60 ± 0.05 



Element 

Mass 

Number 

[X] 

fraction 



H 

9.99 x 10" 1 


1.0 

0.132 

He 

<1.98 x 10“ 5 

<5.0 x 

10“ 6 

<-4.099 

C 

7.15 x 10“ 6 

6.0 x 

10- 7 

-2.520 

N 

2.08 x 10“ 6 

1.5 x 

10- 7 

-2.522 

O 

1.90 x 10“ 5 

1.2 x 

10“ 6 

-2.479 

Al 

1.12 x 10“ 5 

4.2 x 

10- 7 

-0.695 

Si 

5.01 x 10“ 5 

1.8 x 

10“ 6 

-1.123 

P 

1.54 x 10“ 6 

5.0 x 

10“ 8 

-0.579 

S 

5.72 x 10“ 6 

1.8 x 

10- 7 

-1.733 

Ti 

<3.13 x 10“ 7 

<1.1 X 

10- 7 

<-0.100 

Cr 

<1.66 x 10“ 6 

<5.1 x 

10- 7 

<-0.100 

Mn 

<1.08 x 10“ 6 

<3.2 x 

10- 7 

<-0.100 

Fe 

4.98 x 10“ 4 

9.0 x 

10“ 6 

-0.414 

Co 

<4.92 x 10 -7 

<1.2 x 

10- 7 

<-0.100 

Ni 

3.49 x 10“ 5 

6.0 x 

10- 7 

-0.310 

Zn 

1.30 x 10“ 5 

2.0 x 

10- 7 

0.873 

Ge 

3.24 x 10“ 6 

4.5 x 

10“ 8 

1.135 

Sn 

3.53 x 10“ 7 

3.0 x 

10“ 9 

1.569 

Eb-v 

0.0005 ± 0.0005 




log (V H , /cm 2 ) 

18-34!“ 8 




log (AW cm 2 ) 

13.54!“;“ 




4/ pc 

62 + 4 




M / M 0 

0 555 +0035 

U.J^_0 029 




R/Rq 

0 0195 +0 0004 

U,U17J - 0 . 00 05 




log (L / L B ) 

0 63 +0 ‘ 37 

-0.34 





the Galactic plane of z = 8 ± 1 pc. This is in agreement with 
the HlPPARCOS 17 parallax measurement (van Leeuwen 2007, 
HIP23692) of d — 59.88_^° 3 pc and the XHIP 18 value of 
d = 57.96 ± 10.31 pc (Anderson & Francis 2012). The spec- 
troscopic distance of d = 55.84± 0.86pc determined by Flolberg 
et al. (2008) is slightly smaller, this error estimate, however, ap- 
pears to be too optimistic. 


17 http : //www.rssd. esa. int/index .php?project=HIPPARCOS 

18 Extended HlPPARCOS compilation 
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Fig. 26. Comparison of our final TMAP model (blue, thick) with a 
TMAW model (red, thin). Top : astrophysical fluxes at the stellar sur- 
face, middle', ratio TMAP/TMAW flux, bottom: normalized fluxes. 


5. TheoSSA: synthetic stellar spectra on demand 

At the end, we want to compare our final TMAP model flux with 
an SED, that was calculated with the TMAW tool (Sect. 1) and 
considers H, He, C, N, and O only. Figure 26 shows that the 
TMAP flux is higher everywhere at A > 911 A. The reason is 
strong metal-line blanketing at A < 91 1 A that causes a flux in- 
crease at longer wavelengths. It amounts to about 10% at 1000 A 
and to about 5% at 7000 A. The lower panel of Fig. 26 illustrates 
that the theoretical line profiles of the Hi Balrner series are al- 
most identical in both, TMAP and TMAW models, with the ex- 
ception of an increased emission reversal in the line core of Her 
in the TMAP model. 

The TheoSSA database contains currently TMAP SEDs of 
a dozen standard stars. Some of them are represented by differ- 
ent models for comparison, because they were initially calcu- 
lated for the calibration (Vernet et al. 2008a,b, 2010) of ESO’s 19 
second generation VLT 20 instrument XSHOOTER 21 (Vernet 
et al. 2011) while the parameters of Gianninas et al. (201 1) and 
Giammichele et al. (2012) were published later (Table 9). 


19 European Southern Observatory. 

20 Very Large Telescope. 

21 http://www.eso.org/sci/facilities/paranal/ 
instrument s/xshooter 


Table 9. Standard star SEDs (references for T e ft and logg are given) 
presently available in TheoSSA. 


Name 

WD no. 

Spectral type 

Teff 

[K] 

log g 

[cm/s 2 l 

EG 274 

1620-391 

DA2 (+ G5V) 

24276 

8.01" 




25 980 

7.96*" 

Feige 67 


Op+WDsd 

75 000 

5.20' 1 

Feige 110 

2317-054 

sdO 

40000 

5.00" 

G191-B2B 

0501+527 

DAO 

58 883 

7.46" 




61 193 

7.49* 




60920 

7.55* 




60000 

7.55» 

G 93-48 

2149+021 

DAZ3 

18100 

7.85* 




18170 

8.01* 

GD50 

0346-011 

DA2 

40550 

9.22' 




42 700 

9.20* 

GD71 

0549+158 

DAI 

32 747 

7.68* 




32 780 

7.83 ; 




33 590 

7.93* 

GD 108 

0958-073 

sdB 

22908 

5.30* 

GD 153 

1254+223 

DAI. 5 

38 205 

7.89' 




38 686 

7.66* 




40590 

7.93* 

HZ 2 

0410+117 

DA3 

20600 

7.90* 




21600 

7.98* 

HZ43A 

1314+293 

DAl+dM3e 

51 116 

7.90* 




56 800 

7.89* 

Sirius B 

0642-166 

DA2 

24 826 

8.60* 




25 970 

8.57*-" 


Notes. WD numbers are from McCook & Sion (1999). (o) Assumed; 
(fc) Gianninas et al. (2011); (c) Giammichele et al. (2012); (d> Bauer & 
Husfeld (1995); (e) Heber et al. (1984a), He mass fraction of 0.107; 
O Finley et al. (1997); (9) this work; w Guseinov et al. (1983); 
® Barstow et al. (2001 ); (J> Kilkenny et al. (1988); (i) Beuermann et al. 
(2006). 


6. Accuracy of flux calibration with G191-B2B 

Our spectral analysis was performed using state-of-the-art 
atomic data and model-atmosphere code. The best reproduction 
of UV and optical spectra was achieved with chemically homo- 
geneous models. In Fig. 27, we compared two models at the edge 
of our error ranges in T e ff and log g. The deviation in the con- 
tinuum flux of two TMAP model SEDs is ~3% in the optical 
and ~5% in the FUV. A systematic error is present due to the 
uncertainty of the used atomic data, such as oscillator strengths 
where it is typically ~15% for a single line. The employment of 
many lines of many ions of many atoms in a spectral analysis 
minimizes the propagation of these uncertainties into the errors 
of the main photospheric properties like T ei f, log g, and the abun- 
dances. An additional systematic error may be present between 
individual model-atmosphere packages (e.g. Rauch 2008b) be- 
cause of differences in coding, approximations, etc. 

The situation for the flux calibration is, however, not that 
strongly dependent on the exact Teff and log g values (the latter 
is even less important). E.g. in the case of G191-B2B and our 
errors (3% in T e a and 0.05 dex in logg), a normalization to 
a precisely measured brightness will reduce the deviation be- 
tween model SED and observation much below 1% in the op- 
tical and infrared. The residuals among the three primary stars 
G191-B2B, GD71, and GD 153 are generally sub percent at the 
longer wavelengths (Bohlin 2007). The remaining deviation in 
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Fig. 27. Flux ratio of two TMAP models at the edge of the r eff and 
log g error ranges and our final model. 



year 

Fig. 28. Determined T e g and log g values in the last 34 years. Blue tri- 
angles denote analyses with LTE models, red squares those with NLTE 
models (Table 2). The result of Koester et al. (1979, log g = 5.95) is out- 
side the top and bottom panels. 


large scatter in the 1990s (error ranges are not shown for clar- 
ity), the three most recent analyses, that are all based on so- 
phisticated NLTE model-atmosphere techniques, show a good 
agreement within relatively narrow error ranges of about 3% in 
T e e and 0.05 dex (~12%) in log g. Ironically, these latest results 
agree quite well with the very first line-profile analysis presented 
by Holberg et al. (1986) performed with a Lycr line-profile fit 
with a simple, pure-H LTE model atmosphere. The TheoSSA 
database may help to get closer to the intended goal of 1% accu- 
racy in absolute flux calibration. 

We presented here our spectral analysis of G191-B2B to 
demonstrate the current state-of-the-art. We are presently able 
to reproduce the observed spectrum from 250 A to the infrared. 
The EUV part from 150 A to 250 A cannot be modeled, neither 
by our homogeneous nor by our stratified models. The reason is 
unknown. 

A similar analysis of the UV spectrum of the calibration star 
BD+28°4211 (T e ff = 82 000 ± 5000 K, \ogg = b.2 + _° 0 \) was just 
published by Latour et al. (2013). 

Model-atmosphere codes have arrived now at a high level of 
sophistication, and we already encounter problems getting re- 
liable atomic data to reproduce the high-resolution and high- 
S/N spectra that are obtained with presently available instru- 
ments. This is a challenge for atomic physicists to be prepared 
for upcoming telescopes and instruments. 
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the FUV wavelength range is presently less than 2%. Further 
improvement is essentially dependent on the reliability of the 
atomic data (Sect. 7). 

7. Conclusions 

The TheoSSA service is designed to provide theoretical stel- 
lar SEDs of any kind in VO-compliant form. Its efficiency is 
strongly increasing if more different model-atmosphere groups 
provide their SEDs with a proper description in their respec- 
tive meta data. The establishment of a database of spectropho- 
tometric standard stars is an opportunity to use the same 
model SEDs for astrophysical flux calibration. Many model- 
atmosphere groups have their own best models for some of 
these stars, for which a common base for comparison arises. 
Differences in the algorithms for considered physics, assump- 
tions, and approximations in different model-atmosphere codes, 
lead to systematic deviations in general. 

Figure 28 shows the temporal development of the T e g and 
log g determinations of G191-B2B. While both values had a 
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Fig. 1 . Flow diagram of TheoSSA. The VO user sends an SED request to the GAVO database by entering the photospheric parameters. If a suitable 
model is available within the desired tolerance limits, it is offered as a results table. In case that the parameters are not exactly matched, the VO user 
may decide to calculate a model with the exact parameters. TMAW will start a model-atmosphere calculation at our institute’s (IAAT) PC cluster 
then. Extended model grids make use of computer resources that are provided by AstroGrid-D. As soon as the model is converged, the VO user 
can retrieve the SED table from the GAVO database. 
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IV . ... .. . VIII log ( m / g/cm 2 ) r eff = 60000K log (g / cm/s 2 ) = 7.60 


Fig. 6. Temperature and density structure and ionizations fractions of our model with T e g = 60 000 K and log tj = 7.60. IG denotes a generic model 
atom consisting of Ca, Sc, Ti, V, Cr, Mn, and Co. 
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Fig. 10. Synthetic line profiles of C III- IV, N III- V, Si III- IV, PlV- V, S IV- V, and GelV-V calculated from models with 57 000 K < T^g < 62 000 K and log g = 7.6 compared with FUSE and STIS 
observations. The small horizontal bars, labeled with the element’s name, indicate T t g where the ionization balance is best reproduced. The lines are identified at the top. “is” denotes interstellar 
origin. In the cases of C IV dd 1548, 1550 A, C IV AA 1548, 1550 A, and Si III AA 1206.500, 1206.555 A, the photospheric line profile is shown by a dashed, blue line (for the T c g = 57 000 K model at 
the top only). 
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Fig. 11. Same as Fig. 10 for FelV- VI lines (from left to right panels, marked blue in the top panels) only. 
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Fig. 12. Same as Fig. 10 for NilV- VI lines (from left to right panels, marked green in the top panels) only. 
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Fig. 23. Comparison of our model spectra (all calculated with logg = 7.60) with the STIS observation. Top: final homogeneous model (TMAP). 
The stratified NGRT models have = 65 000 K, 60 000 K, 55 000 K, and 50 000 K (from top to bottom). The lines in the models are marked. 
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Table 2. T e g and log g from previous analyses of G191-B2B. 


T e ff 

log g 

Reference 

Method 

[K] 

[cm/s 2 ] 



61900 

7.5" 

Shipman (1979) 

LTE, pure H, optical colors 

56 788 ± 3336 

5.95 + 0.04 

Koester et al. (1979) 

LTE, pure H. optical colors,* 

62 250 ± 3520 

7.55 ± 0.35 

Holberg et al. (1986) 

LTE, pure H, H I L a line 

59 250 ± 2000 

7.50 ±0.10 

Kidder (1990) 

cited by Holberg et al. (1991) 

LTE, pure H. H I Hy and H<S lines 

61000!““ 

8.00" 

Green et al. (1990) 

LTE, H+He, extreme ultraviolet (EUV) continuum 

61170*®! 

8.00" 

Finley et al. (1990) 

LTE, H+He, ultraviolet (UV) continuum 

62 250+ 1000 

8.00" 

Finley et al. (1990) 

LTE, H+He, Hi Lcr line 

62 250 

7.55 

Vennes et al. (1991) 

NLTE, H+HeCNSi, UV spectrum 

53 500 ± 500 


Koester & Finley (1992) 

LTE, H+He" , UV continuum 

60500 + 900 

7.50 ± 0.05 

Vidal-Madjar et al. (1994) 

LTE, H+He, Hi Balmer lines 

57 900+ 1500 

7.50" 

Dupuis et al. (1995) 

NLTE, pure H, EUV continuum' 2 

54000 + 800 

7.50" 

Dupuis et al. (1995) 

NLTE, H+CNOFe, EUV continuum'' 

60500+ 1000 

7.5 

Lanz et al. (1996) 

LTE, H+HeC, Hi Balmer lines 

56 000 ± 1000 

7.5 

Lanz et al. (1996) 

NLTE, H+HeC, H I Balmer lines 

55 200 ± 1000 

7.5 

Lanz et al. (1996) 

NLTE, H+HeCFe, Hi Balmer lines 

64000+ 1000 

7.64 ± 0.06 

Vennes et al. (1996) 

LTE, pure H, H 1 Lyman lines 

57 900+ 1500 

7.5 

Vennes et al. (1996) 

LTE, pure H, extreme UV (EUV) continuum* 

64100 + 700 

7.69 ± 0.04 

Vennes et al. (1996) 

LTE, H+CNOFe, H 1 Balmer lines 

52 600 ± 800 

7.53 ± 0.07 

Vennes et al. (1996) 

LTE, H+CNOFe, Hi Lyman lines 

54000 + 800 

7.5 

Vennes et al. (1996) 

LTE, H+CNOFe, EUV continuum* 

61 193 + 241 

7.49 ± 0.01 

Finley et al. (1997) 

LTE, H-Ni, H I Balmer lines 

59 160 +1270 

jy auu_ 10 70 

7 36 +0 08 

' " J -0.07 

Barstow et al. (1998) 

NLTE, pure H, Hi Balmer lines 

59 190!““ 

7 36 +0 - 07 

' • JU -0.07 

Barstow et al. (1998) 

NLTE, H+He, H 1 Balmer lines 

59060 +113 ° 

jy uuu _ 1090 

7 36 +0 08 

' ' J -0.07 

Barstow et al. (1998) 

NLTE, H+He + heavy-metal poor. Hi Balmer lines 

53 840 ! 9 “ 

7 QQ+0-07 
/,JO -0.08 

Barstow et al. (1998) 

NLTE, H+He + heavy-metal rich, H 1 Balmer lines 

52 920 ± 350 

7.36 ± 0.03 

Barstow et al. (1998) 

NLTE, H+He + heavy-metal rich. Hi Lyman lines" 

56000 

7.6 

Wollf et al. (1998) 

LTE + NLTE, H+CNOSiFeNi, EUV continuum* 

56000 

7.6 

Dreizler & Wolff (1999) 

NLTE, H+CNOSiFeNi, diffusion model, EUV to optical 

54600 + 200 

7.60 ± 0.02 

Barstow et al. (2001) 

NLTE, H+HeCNOSiFeNi, Hi Balmer lines 

52930 + 3600 

7.16 ± 0.2 

Barstow et al. (2001) 

NLTE, H+HeCNOSiFeNi, Hi Lyman lines' 

53 180 ±530 

7.43 ± 0.04 

Barstow et al. (2001) 

NLTE, H+HeCNOSiFeNi, Hi Lyman lines" 

56000 

7.59 

Schuh et al. (2002) 

NLTE, H+HeCNOSiFeNi, diffusion model, EUV continuum* 

54000 

7.5 

Holberg et al. (2003) 

NLTE, metal lines 

58 865 ± 706 

7.57 ± 0.038 

Lajoie & Bergeron (2007) 

NLTE 9 , pure H, H I Balmer lines 

60680+ 15 000 

7.57* 

Lajoie & Bergeron (2007) 

NLTE 9 , pure H, H I Lyman lines' 

57 414 + 4700 

7.57* 

Lajoie & Bergeron (2007) 

NLTE 9 , pure H. V-normalization method 

61980 + 514 

7.56 ± 0.04 

Allende Prieto et al. (2009) 

NLTE, H' , H I Balmer lines 

60920 + 993 

7.55 ± 0.05 

Gianninas et al. (2011) 

NLTE, H+CNO, H I Balmer lines* 


Notes. Assumed log g value. (i) The authors note that the results are extrapolated from their model grid. w Stratified model, hydrogen- layer mass 
between 6 x 10 -15 and 8 x 10 -15 M 0 . ( * ) Extreme Ultraviolet Explorer (EUVE, http://heasarc.gsfc.nasa.gov/docs/euve/euve.html) 
observations. Orbiting and Retrievable Far and Extreme Ultraviolet Spectrometer (ORFEUS, http://www.uni-tuebingen.de/en/4221) 
and FUSE observations. ^ Hopkins Ultraviolet Telescope (HUT, http://praxis.pha.jhu.edu/) observations. (9) International Ultraviolet 
Explorer (IUE) observations. ( *’ Adopted from their optical solution. w Models described in Liebert et al. (2005). Ij> The authors note that T e g 
may be overestimated by »6000 K because their pure-H models are inappropriate due to the photospheric metal content. ® New H I Stark line- 
broadening tables from Tremblay & Bergeron (2009). 
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Table 3. Abundances of photospheric trace elements in G191-B2B 
from previous analyses. 


Element 

Log mass fraction 

He 

-4.2 ±0.1° 
-4.4 ± 0.3* 

-4.4' 

C 

-4.6 ± 0.3" 
-4.6 ± 0.3* 

-5.6' 

N 

-4.3 ± 0.4'' 
-5.6 ± 0.3* 

-5.9' 

O 

-4.8 ± 0.3'' 
-4.8 ± 0.3* 

-4.6 ' 

Mg 


<— 5.6' 

Al 


-5.1' 

Si 

-5.1 ± 0.4'' 
-5.1 ±0.5 C 
-5.0 ± 0.3* 

-5.0' 

P 

-6.2 ± 0.2' 


S 

-5.2 ± 0.5 C 


Cl 

<-7.0* 


Cr 


<— 6.3' 

Mn 


<— 6.3' 

Fe 

-3.8 ±0.3* 
-3.4 ± 0.4-* 
-3.3 ± 0.3* 

-3.8' 

Co 


<— 6.2' 

Ni 

-4.2 ± 0.5 9 
-4.2 ± 0.4^ 
-3.9 ±0.3* 

-4.4' 

Ge 

-6.1 ± 0.2* 


Sn 

-6.9 ± 0.2* 



Notes. He abundance assumed, no abundance uncertainties given. 
(a) Cruddace et al. (2002); (i>) Lanz et al. (1996); Vennes et al. (1996); 
w Vidal-Madjar et al. (1994); (e) Chayer et al. (1996); Werner & 
Dreizler (1994); <9) Holberg et al. (1994); m Vennes et al. (2005); 
(l) Holberg et al. (2003). 
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Table 4. Statistics of our model atoms. 


Ion 

NLTE levels 

LTE levels 

Lines 

Sample lines 

Hi 

14 

2 

91 


II 

1 

- 

- 


He I 

29 

74 

69 


II 

16 

16 

120 


III 

1 

- 

- 


Cn 

1 

45 

0 


III 

44 

23 

190 


IV 

54 

4 

295 


V 

1 

0 

0 


N II 

1 

246 

0 


III 

34 

32 

129 


IV 

90 

4 

546 


V 

54 

8 

297 


VI 

1 

0 

0 


On 

1 

46 

0 


III 

72 

0 

322 


IV 

38 

56 

173 


V 

76 

50 

472 


VI 

54 

8 

291 


VII 

1 

0 

0 


A1 II 

1 

4 

0 


III 

7 

29 

10 


IV 

6 

183 

3 


V 

6 

223 

4 


VI 

1 

0 

0 


Sim 

17 

17 

28 


IV 

16 

7 

44 


V 

1 

0 

0 


Pill 

3 

7 

0 


IV 

21 

30 

9 


V 

18 

7 

12 


VI 

1 

0 

0 


Sill 

1 

230 

0 


IV 

17 

83 

32 


V 

39 

71 

107 


VI 

25 

12 

48 


VII 

1 

0 

0 


Fe III 

7 

0 

25 

537689 

IV 

7 

0 

25 

3 102 371 

V 

7 

0 

25 

3 266247 

VI 

7 

0 

33 

991 935 

VII 

7 

0 

39 

200455 

VIII 

1 

0 

0 

0 

Ni III 

7 

0 

22 

1 033 920 

IV 

7 

0 

25 

2512561 

V 

7 

0 

27 

2766664 

VI 

7 

0 

27 

7 408 657 

VII 

7 

0 

33 

4 195 381 

VIII 

1 

0 

0 

0 

IG ill 

1 

0 

0 

0 

IV 

7 

0 

25 

1579918 

V 

7 

0 

23 

2230921 

VI 

7 

0 

25 

1455 521 

VII 

7 

0 

24 

1 129512 

VIII 

1 

0 

0 

0 


Notes. IG denotes a generic model atom consisting of Ca, Sc, Ti, V, Cr, Mn, and Co. “Sample lines” are individual Kurucz lines that are sampled 
to superlines for Ca-Ni (Rauch & Deetjen 2003). 
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Table 4. continued. 


Ion 

NLTE levels 

LTE levels 

Lines 

Sample lines 

Zn II 

1 

5 

0 


III 

2 

10 

0 


IV 

31 

0 

87 


V 

5 

15 

0 


VI 

1 

0 

0 


Ge III 

1 

15 

0 


IV 

8 

1 

8 


V 

85 

0 

878 


VI 

11 

25 

0 


VII 

1 

0 

0 


Sn III 

3 

18 

2 


IV 

6 

4 

1 


V 

5 

4 

0 


VI 

6 

0 

0 


VII 

1 

0 



Total 70 

1038 

1614 

4646 

32411752 
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Fig. A.l. Comparison of the FUSE observation with our final model. Stellar and interstellar lines are identified, “i 
“unid.” denotes unidentified lines. 



i” denotes interstellar origin, 
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Fig. A.2. Comparison of a section of the STIS observation with our final model. Stellar and interstellar lines are identified, “is” denotes interstellar 
origin, “unid.” denotes unidentified lines. 
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